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Abstract: Finite volume renormalization scheme is one of the most fascinating scheme 
for non-perturbative renormalization on lattice. By using the step scaling function one can 
follow running of renormalized quantities with reasonable cost. It has been established the 
Schrodinger functional is very convenient to define a field theory in a finite volume for the 



Q_i! renormalization scheme. The Schrodinger functional, which is characterized by a Dirichlet 

boundary condition in temporal direction, is well defined and works well for the Yang-Mills 



theory and QCD with the Wilson fermion. However one easily runs into difficulties if one 
sets the same sort of the Dirichlet boundary condition for the overlap Dirac operator or 
the domain-wall fermion. In this paper we propose an orbifolding projection procedure 
to impose the Schrodinger functional Dirichlet boundary condition on the domain-wall 
fermion. 



1. Introduction 

Perturbative renormalization factor is a source of systematic errors in numerical investi- 
gation of lattice QCD. There has been progress in numerical simulation with dynamical 
fermions nowadays and sources of systematic error is decreasing. Evaluation of renormal- 
ization factors in non-perturbative method is also required. Finite volume renormalization 
scheme is one of the most fascinating procedure to define non-perturbative renormalization 
scheme on lattice. By using the step scaling function one can follow running of renormalized 
quantities from low energy region to perturbative region with reasonable cost for recent 
computers. It has been established that the Schrodinger functional is very convenient to 
define a field theory in a finite volume for renormalization scheme. 

The Schrodinger functional (SF) is defined as a transition amplitude between two 
boundary states with finite time separation [1, 2, 3, 4] 

Z = (C';x = T\C;x = 0) = f V^e~ S[ ^ (1.1) 



and is written in a path integral representation of the field theory with some boundary 
condition. For renormalization of a finite volume theory defined through the SF the renor- 
malization scale is introduced by a finite volume TxL 3 ~ L 4 of the system. The formulation 
is already accomplished for the non-linear a- model [5], the non-Abelian gauge theory [6] 
and the QCD with the Wilson fermion [7, 8] including 0{a) improvement procedure [9, 10]. 
(See Ref. [11] for review.) 

Several renormalization quantities like running gauge coupling [12, 13, 14, 15, 16, 
17, 18], Z-factors and 0{a) improvement factors [19, 20, 21, 22, 23, 24] are extracted 
conveniently by using a Dirichlet boundary conditions for spatial component of the gauge 
field 

A k (x)\ X0=0 = C fc (f), A k (x)\ X0=T = C' k (x) (1.2) 

and for the quark fields 

P+ip(x)\ xo=0 = p(x), P-ip(x)\ X0=T = p'(x), (1.3) 

ip(x)P-\ X0=0 =p(x), 4>(x)P + \ X0= t = p'(x), (1.4) 

ft = ^p- (1.5) 

One of advantage of this Dirichlet boundary condition is that the system acquire a mass 
gap proportional to 1/T and there is no infra-red divergence. The finite volume plays a 
role of an infra-red cut-off. Field theory with Dirichlet boundary condition is shown to be 
renormalizable for the pure gauge theory [6] and QCD with the Wilson fermion [8]. 

Although it is essential to adopt Dirichlet boundary condition for a mass gap and 
renormalizability, it has a potential problem of zero mode in fermion system. For instance 
starting from a free Lagrangian 

£ = V (Tm^m + m ) ^ (!-6) 



with positive mass m > and the Dirichlet boundary condition 

P-V>|xo=0 = 0, P+VUo=T = (1.7) 

the zero eigenvalue equation (70^0 + m ) V' = i n temporal direction allows a solution 

V> = P+e- mX(> + p_ e - m ( T -*») (1.8) 

in T — > 00 limit and a similar solution remains even for finite T with an exponentially 
small eigenvalue oc e~ . In the SF formalism this solution is forbidden by adopting an 
"opposite" Dirichlet boundary condition (1.3) and the system has a finite gap even for 
m = [7]. 

For the Wilson fermion [7] on lattice the Dirichlet boundary condition is automatically 
chosen among 

P±^|x =0 = 0, P T Mx =T = (1.9) 

depending on signature of the Wilson term. For example if we adopt a typical signature of 
the Wilson term 

Dw = 1»\ (K + V m) " f V*V„ + M ( L1 °) 

the allowed Dirichlet boundary condition is the same as (1.3). In this case the zero mode 
solution is forbidden by choosing a proper signature for the mass term; the mass should be 
kept positive M > to eliminate the zero mode [7]. 

However as was discussed in the previous paper [25, 26] this zero mode problem may 
become fatal in the overlap Dirac operator [27, 28] and the domain- wall fermion [29, 30, 31]. 
Both the overlap Dirac operator and the domain-wall fermion is defined through the four 
dimensional Wilson Dirac operator (1.10) but with an opposite signature for the Wilson 
fermion mass parameter M (domain-wall height) to the Wilson parameter r. An opposite 
signature is necessary to impose heavy masses on the doublers and a single massless mode 
to survive. A requirement to the four dimensional Wilson Dirac operator is that Dyy should 
not have a continuous zero mode. If this is not the case the chiral Ward-Takahashi identity 
is broken dynamically for the domain-wall fermion that the explicit breaking term does not 
vanish [31]. For the overlap Dirac operator a continuous zero mode may break locality of 
the Dirac operator [32]. 

If the Dirichlet boundary condition (1.3) (1.4) is imposed to all fermion fields of the 
overlap Dirac operator or the domain-wall fermion exponentially small eigenvalues are 
allowed in the kernel Dw because of an opposite signature of the Wilson parameter and 
the domain- wall height. Since these small eigenvalues are continuous in spatial momentum 
they may be a lethal problem in large T limit to break essential properties of the chiral 
Dirac operator. 

One may wonder that the small eigenvalues are boundary effect and should be localized 
near the temporal boundary. If one considers physics apart from the boundary there should 
be no harm. However this is not the case for our purpose to define renormalization scheme. 



In finite volume scheme the renormalization scale is given by a size of the box, which is 
realized by considering a correlation function of operators to be separated by an order of 
box size. At least one of operators cannot be away from the boundary. Furthermore it is 
convenient for the SF scheme to set one of the operator at the boundary. 

In order to solve this problem an orbifolding projection procedure was proposed for 
the overlap Dirac operator in Ref. [25]. * In this formulation we start from a theory on 
S 1 x R 3 and impose orbifolding projection S 1 jZ% on temporal direction. Since we have 
set anti-periodic boundary condition in temporal direction S 1 before projection we have 
a mass gap proportional to 1/T, which is not broken by the orbifolding. Because of this 
mass gap we can avoid the zero mode problem of Dirichlet boundary condition. 

In this paper the orbifolding formulation of the SF boundary condition is applied to 
the domain-wall fermion. In section 2 the domain-wall fermion on S* 1 x T 3 is introduced. 
Formulation of domain-wall fermion in finite volume with the SF boundary condition is 
discussed in section 3. Application of orbifolding procedure to fermionic part is almost 
straightforward as was discussed in Ref. [26]. We can use the same kind of symmetry 
argument as in the previous paper [25]. Difficulty is in a treatment of the Pauli-Villars 
field. We adopted effective Dirac operator for this purpose. The proper Dirichlet boundary 
condition (1.3) (1.4) may not be the unique choice to define a finite volume renormalization 
scheme. In section 4 a chirally twisted boundary condition is discussed to define a finite 
volume field theory keeping a good property of the SF boundary condition. Section 5 is 
devoted for conclusion. 

2. Domain- wall fermion action 

The purpose of this paper is to introduce the domain-wall fermion system, with which 
we can define a finite volume renormalization scheme (Schrodinger functional scheme). 
The formulation for the pure Yang-Mills theory has been established in Ref. [6] by using 
a transition amplitude between two boundary states (Schrodinger functional). In this 
formulation the gauge field (link variable) lives in a finite box Nt x N^ with a periodic 
boundary condition in spatial direction and the SF Dirichlet boundary condition at the 
temporal boundary 

U k (x,0) = W k (x), U k {x,N T ) = W' k (x). (2.1) 

We shall adopt this procedure for the gauge part and treat the gauge field as an external 
field in this paper. 

The transition amplitude of the fermion field has been introduced for the Wilson 
fermion using the transfer matrix in Ref. [7]. The fermion field resides in the same finite 
box for the path integral formalism with periodic or twisted boundary condition [9] in 
spatial direction and the SF Dirichlet boundary condition (1.3) and (1.4) in temporal 
direction. This fermion system is renormalizable including a shift in the boundary field p 



'After finishing this paper a new paper appeared to propose a method to define chiral symmetric theory 
with the SF Dirichlet boundary condition [33]. 



and p [8]. Another specific property is that this system has a mass gap proportional to the 
temporal length 1/T and the finite box serves as an infra-red regulator. 

We shall construct the domain-wall fermion system in a finite box keeping the same 
sort of properties as the Wilson fermion; (i) the theory has a mass gap proportional to 1/T, 
(ii) there are boundary fields p and ~p in temporal direction and the theory is renormalizable 
including a shift in these fields. If one naively impose the boundary condition (1.3) and (1.4) 
to all the fifth dimensional field ip(x, s) then the chiral symmetry is broken "dynamically" 
as explained in the introduction. In order to avoid this problem we adopt an orbifolding 
procedure, where we start from doubled time length 2Nt and fermion fields in the finite box 
of length Nt with the Dirichlet boundary condition is realized by an orbifolding projection. 
For this purpose we copy gauge configuration with the SF boundary condition (2.1) into 
negative region and produce a time reflection symmetric configuration, which satisfies 

U k (x,x ) = U k (x,-x ), U (x,x ) = U^(x,-x - 1) (2.2) 

as in the previous formulation of overlap Dirac operator [25]. The periodic boundary 
condition is set with length 2Nt 

Up&xo + 2N T ) = U^x,x ). (2.3) 

In this paper we adopt the Shamir's domain-wall fermion [30, 31] on a lattice 2Nt x 
N 3 L xN 5 

N T N 5 

S = Y1 12 Y,^(x,s)D dw{ (x,y;s,t)^(y,t). (2.4) 

x,y xo,Ho=—Nt+1 s,t=l 

xq and j/o represent the temporal coordinate which runs —Nt + 1 < ^o < Nt- s and t 
are used for the fifth dimensional coordinate which runs 1 < s < N$. For later use of 
orbifolding we set the anti-periodic boundary condition in temporal direction 

ip(x,x + 2N T ,s) = -tp(x,x ,s), ip(x,x + 2N T ,s) = -ip(x,x ,s). (2.5) 

The Dirac operator is given as a five dimensional Wilson's one with conventional Wilson 
parameter r = 1 and negative mass parameter (domain-wall height) — M with < M < 2 

Ddwd x > V, s > t) = jmD m - -D 2 - M 

= ( zl ip c/ °(*)H^ + zi y^ c/ o t (y)^o) s Xi , m s s , t 

+ f g — % ~ U *( x ) S Vi>*i+ 1 + 2 — lU i (y^yi'Xi- 1 ) S xo,yo S s,t 

+ ( 5 — -^ + ( m f)s,t H — -£l~(m f ) s j J 5 XjV 

+ (5 - M)S x>y S a>u (2.6) 



where W are hopping operator in temporal direction with anti-periodic boundary condi- 
tion, whose explicit form for 2Nt = 6 is written as 



W. 
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(2.7) 



Q are hopping operator in fifth direction with Dirichlet boundary condition (for massless 
case), whose matrix form for N5 = 6 is given by 



tt + (mf) S! t 



( 1 0\ 

10 

10 

10 

1 

\-m f 0/ 



n-( m/ ) = (n + (m f )y . 



(2.8) 



Here rrif is a physical quark mass. 

The physical quark field is defined by the fifth dimensional boundary field with chiral 
projection 



(2.9) 
(2.10) 

(2.11) 



q(x) ■- 


= (Pl$s,i + PrS s ,n 5 )iP(x, s) 


q(x) -- 


= ip(x, s) (5 StN5 P L + 6 S;1 P R ) 


Po, r 


__ 1±75 



R/L 



The physical quark mass term is given as an ordinary form C r 
field. 



rrifqq with this quark 



3. Schrodinger functional with conventional boundary condition 

In this section we shall construct the domain-wall fermion system in finite box, in which 
the conventional SF Dirichlet boundary condition (1.3) and (1.4) is satisfied by the physical 
quark field. This formulation will be done by making use of an orbifolding in temporal 
direction. 



3.1 Orbifolding construction of SF boundary condition 

Since we adopted (anti) periodic boundary condition in temporal direction with period 
2Nt fields live on S 1 . The orbifolding S 1 /Z2 is to identify the negative time coordinate 
with the positive one xq = —xq. Identification of fields on S 1 is performed according to the 
symmetry of the theory including the time reflection. A homogeneous Dirichlet boundary 
condition will appear at fixed points. 



The time reversal symmetry of the domain- wall fermion is given by 

i){x, x , s) ->• T, XOiyo . Sjt tp(x, y , t), ip(x, xq, s) ->• ip(x, y , t)T, y0iXo;tjS , (3.1) 



^-"xojj/o^jt — ^T5T0-*keojOT *>t' 



(3.2) 



P. 



s,t 



(iV 5 = 6) 



(3.3) 



where P is a parity transformation in fifth direction P Sjt ip(x, xq, t) = tp(x, xo, ./V5 — s + 1), 
whose matrix representation is 

/0 1\ 

10 

10 

10 

10 

\l 0/ 

and i? is a time reflection operator acting on the temporal direction R X(hyo ip(x,yo,s 
ip(x, —xq, s), whose matrix form is given by 

/0 1 \ 

10 

10 

10 

10 

\0 -1/ 

to satisfy anti-periodicity in 2Nt- We notice that R has a symmetric fixed point xq 
and an anti-symmetric fixed point xq = Nt 
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■•''o,;yi) 



(2N T = 6) 



(3.4) 







Rip(x,0,s)=tp(x,0,s), RiI>(x,Nt,8) = -i/>(x,N t ,8). (3.5) 

The domain-wall fermion Dirac operator is invariant under the time reflection 

[S, D dw{ ] = (3.6) 

since the reflection invariant gauge configuration (2.2) is adopted . 

In order to realize the SF boundary condition at the fixed points we need to combine 
the chiral transformation with the time reflection [25]. The chiral transformation is given 
by a vector like rotation of fermion field but with a different charge for two boundaries in 
fifth direction [31] 

ip{x,s) ->■ iQ s ,t^(x,t), ip(x,s) -> -ip(x,t)iQ t ,s, (3.7) 

where Q is the vector charge matrix which flips sign in the middle of the fifth direction 
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(N 5 = 6). 



(3.8) 
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We consider massless rrif = theory in this sub-section. 

Here we should notice that this chiral rotation is not an exact symmetry of the domain- 
wall fermion Dirac operator but we have an explicit breaking term 

QD dwi Q - D dwl = 2X, (3.9) 

where X is a contribution from the middle layer, which picks up a charge difference there 

x = ( p LS s ^S t ^ +1 + P R S s ^ +1 6 t M^j S x , y . (3.10) 

However it was discussed in Ref. [31] that if we consider correlation functions between the 
bilinear ipXip and the physical quark operators contribution is suppressed exponentially 
in N$ under the condition that the transfer matrix in fifth direction has a gap from unity. 
Furthermore the domain-wall fermion Dirac operator with explicit time reflection invariance 
(3.6) does not have index [25], since the contribution to the index [34] 

lim a 4 Y^(^(x,s) l5 X s ^(x,t)) = - lim tr (75X-M (3.11) 

can be shown to vanish by using anti-commutativity {75A, S} = 0. We expect that X has 
no effect on anomaly. We shall ignore this term in the following by constraining that we 
treat the physical quark Green's functions only. 

Another way to avoid the explicit breaking term is to include it into the Dirac operator. 
By using an anti-commutative nature {Q, X} = we can define a chiral symmetric Dirac 
operator by 

D^ = D dw{ + X, (3.12) 

which commutes with Q exactly even at finite N$. The orbifolding projection in the 
following can be defined in an exact sense. A compensation of the exact chiral symmetry 
at finite N§ is a non-locality in the effective Dirac operator, which however is suppressed 
exponentially in N§. Detailed property of this Dirac operator is deferred in appendix A. 

Combining the time reversal transformation (3.1) and the chiral transformation (3.7) 
we define the orbifolding transformation 

ip(x, xq, s) -► A XOjyo . Sjt tp(x, 2/0, t), ip(x, x , s) -> tp(x, y , t)A y0:Xo;t:S , (3.13) 

A Xo ,y ; S ,t = r yOlb(PQ) S ,tRx ,yo- ( 3 - 14 ) 

The domain-wall fermion Dirac operator has time reversal symmetry (3.6) and we assume 
that the chiral transformation is an exact symmetry of the Dirac operator 

[Q,D d wf]=0 (3.15) 

by ignoring effect of the explicit breaking term X or by adopting the symmetric Dirac 
operator. The orbifolding transformation becomes symmetry of the Dirac operator 

[AAiwf] = 0. (3.16) 



In order to show this we may use a relation {P, Q} = 0. 

The operator A satisfies a property A 2 = 1 and can be used to define a projection 
operator. The orbifolding identification of the fermion field is given by projecting out the 
following symmetric sub-space 

n_v(x, s ) = o, (^n_)(s, a ) = o, n ± = ^— . (3.17) 

This projection relates fields in negative region to those in the positive ip(x, — xq, s) = 
jQ'y5PQ^(x,xo,s), which means fields in the negative is not independent. As will be 
discussed in appendix B if we consider non-negative region < xq < T, fields in the bulk 
< xq < T is not constrained. Only the boundary fields obey a projection condition 

P-ip(x, 0, s) = 0, P+ip(x,N T ,s) = (3.18) 

@P-) (2,0,8) = 0, @P+)(x,N T ,s) = (3.19) 

with projection operator 

P± = ^^, r = 7o75 PQ. (3.20) 

The orbifolding projection for the physical quark field is given by picking up the bound- 
ary components from the projected fermion field 

(P L 5 sA + P R 5 S>N5 ) (H_) ajt V(M) = U + q(x) = 0, (3.21) 

^{x, t) (n_) M (S S , N5 P L + S 9tl P R ) = g(x)n_ = 0, (3.22) 

n± = i^, r = 7o ^, (3.23) 

which turns out to be the same condition for the continuum theory in Ref. [25]. The proper 
homogeneous SF Dirichlet boundary condition is provided at fixed points xo = 0, Nt for 
the physical quark fields 

P + q(x)\ X0=0 = 0, P-q(x)\ X0=NT = 0, (3.24) 

q(x)P-\ X0=0 = 0, q(x)P + \ X0=Nr = 0. (3.25) 

The massless orbifolded action is given by projection 

s SF = J2 T^dwfV>, ^ = n+Aiwfn + . (3.26) 

We notice the massless SF Dirac operator -D^ f breaks "chiral symmetry" (3.7) explicitly 

by the projection IT + . However the symmetry breaking effect comes from the projection 

(3.18) (3.19) at the boundary. Ordinary chiral Ward-Takahashi identity [31] is satisfied in 

the bulk < xo < Nt where fields are not constrained. 

Since our orbifolded action is given by projecting onto a symmetric sub-space of the 

theory and the orbifolding symmetry is not broken by anomaly ' the renormalizability is 

kept trivially *. 

^According to a similar discussion to that for chiral index (3.11) we can easily show that orbifolding 
matrix A does not have index. 

"'"Boundary source fields are introduced in later sub-section. 



Our original theory on S* 1 has a gap because of the anti-periodic boundary condition. 
This gap is kept intact after orbifolding, which can be confirmed at tree level. We have a 
Hermiticity relation for the SF Dirac operator 

(^dwf) t =75 J P^ f 75 J P (3.27) 

and this Dirac operator connects the same Hilbert sub-space 

Dllf:H--+H-, H_ = {V>|n_^ = 0}. (3.28) 

It is straightforward to solve the eigenvalue problem numerically at tree level. Here we 
omit the detail but we can easily see that the lowest eigenvalue (a gap) converge to ir/2T 
in the continuum limit, which agrees with that of continuum massless theory [7]. 

As will be discussed in appendix B the bulk part of this projected Dirac operator is 
exactly the same as that of the ordinary domain-wall fermion. The physical quark fields 
satisfies the proper boundary condition. Together with the renormalizability and existence 
of the mass gap this orbifolded system is a strong candidate of QCD with the SF boundary 
condition to define a finite volume scheme. 

We have a comment on mass term. We dropped quark mass term since it breaks the 
chiral symmetry. However as was discussed in Ref. [25] it is possible to introduce a mass 
term which is consistent with the orbifolding symmetry (3.13). One of candidates is 

S mass = '^2mfq(x)ri(x )q(x), (3.29) 

X 

where r] is an anti-symmetric step function 

v(-xo) = ~v{xo), r){xo + IT) = 7](x ), 

rt(x ) = 1 for < x < N T . (3.30) 

3.2 Free propagator 

In order to check that the orbifolded domain-wall fermion system describes the QCD with 
the SF boundary condition properly we consider the physical quark propagator at tree level. 
The massless fermion propagator is given as an inverse of the projected Dirac operator 

GZf(x,y,s,t) =2(DZ { )~] L = 2(n + -^!l + ) . (3.31) 



where inverse is defined in the sub-space 7i. 



DZt (DZ)- 1 = (DZtV 1 DZt = n + . (3.32) 
At tree level this propagator can be written in a simple form as 

GZt (x, |/i^) = p£ e^-^GZi (Pi so, Vo; s, t), (3.33) 

L p 

<7 SF (n-Tn ?m- s t) — V^ II J ( p ipo(xo-Vo) i Jpo(x +yo)\ (~p,\ 

G dwi (p,x ,y , S ,t) - 2aNT ^^ DdM ) st/ {{e +e ) [P + ) t , t 

+ Lipo(x -yo) _ e i P o(x +yo)\ (p_^ \ ( 3 _ 34 ) 



where the projection operator P± is defined in (3.20). The temporal momentum po satisfies 
the quantization condition 

2n — 1 
Po = ~^^ n ^ ~ N T + l<n<N T (3.35) 

for anti-periodicity in 2Nt- Ddwf(p) is the domain- wall fermion Dirac operator in momen- 
tum space without orbifolding projection 

aAlwf (p) = iln sin ap^ + W{p) - P L tt + - P R il~ , (3.36) 

W = 1 - M + ^2 (1 - cos ap^) . (3.37) 

The explicit form of its inverse can be derived according to Ref. [30], which we defer to 
appendix C. 

The physical quark propagator is given by selecting the contribution from the boundary 
fields in fifth direction 

GqLk(^ y) = (Pl6 s ,1 + Pr8,,N s ) GJj£f (x, y, S, t) (6 tiNB P L + ^iPfl) 

= 2 (n_G quark n+) :Ci2/ , (3.38) 



where 



G quaTk (x, y) = (P L S aA + PrS s , N5 ) (j^—) (S tyN& P L + S tyl P R ) (3.39) 

V-^dwf / x,y;s,t 

is the physical quark propagator in 2Nt x AT| space-time without any projection. The 
proper Dirichlet boundary conditions [10] 

P+Gll aik (x,y)\ xo=0 = 0, P-G s q ^ rk (x,y)\ Xo=NT =0, (3.40) 

GqLk(^ v)\yo=oP- = 0, Gg^kfo ^Ivo^T^ = (3.41) 

are satisfied for this quark propagator because of the projection H±. By ignoring sub- 
leading terms in e~ Ns the propagator takes the following form at tree level 

x I ( e -iPoJ/o + giPoi/o) p + + ( e -iP0J/0 _ e ipoJ/o) p } ; (3.42) 

which can be shown to approach to the continuum SF propagator of Ref. [10] without any 
0{a) term. This system has no extra zero mode we encountered in naive formulation and 
we conclude that this is equivalent to the QCD with SF boundary condition. 

3.3 Surface term 

In the orbifolding construction of the SF formalism only the homogeneous boundary con- 
dition (3.24) (3.25) can be introduced. However in general SF formalism the Dirichlet 
boundary condition is inhomogeneous as (1.3) and (1.4). The boundary values p,---,p' 



10 



are regarded as external source fields coupled to the dynamical fields and the correlation 
functions involving the boundary fields 

«*> = im- m = -jkr (3 - 43) 

are used conveniently to extract the renormalization factors. 

Coupling between the boundary source fields and the dynamical fields are not intro- 
duced automatically in our formulation since the boundary value vanishes by projection. 
In other words our formulation is protected by an orbifolding symmetry. The action (3.26) 
is invariant under a transformation 

S (n+V) (x, s) = a (TT + V0 (x, s), (3.45) 

S(^U+)(x,s) = -a(^U+) (x,s), (3.46) 

where remaining degrees of freedom il-^ and ip Il_ are kept intact. The boundary source 
fields are elements of U-ip and ip II_. Since the surface term is given to connect the 
boundary source fields and the dynamical fields H+ip, ip n + it is not consistent with the 
orbifolding symmetry. 

In this paper we define a surface term as an orbifolding symmetry breaking term, which 
is consistent with other symmetries of the orbifolded domain-wall fermion. The symmetries 
are; parity 

ip(x,s) -^^ P st ip(-x,xo,t), i/}(x,s)->i/}(-x,xo,t)'YoPts, (3.47) 

charge conjugation 

iP(x,s)^CP st ^ T (x,t), ^(x,s)^ip T (x,t)P ts (-C- 1 ), C = 727 o (3.48) 

and chiral symmetry (3.7) in the bulk < xq < Nt, where chiral Ward-Takahashi identity 
of Ref. [31] is satisfied. 

Using the orbifolding projection (3.17) it is easy to show that the orbifolding transfor- 
mation (3.45) (3.46) is a "chiral" transformation only at the boundary 

5 (P+ip) (x, 0, s) = a (P+tp) (x, 0, s), (3.49) 

5 (^ P+) (x, 0, s) = -a (^ P+) (x, 0, s), (3.50) 

5 (P_V>) (x , N T , s) = a (P_V) (x, N T , s), (3.51) 

S (^ P_) (x, N T , s) = -a (^ P_) (x, N T , s) (3.52) 

and is a vector U(l) transformation in the bulk < xq < Nt 

Sip(x, xq,s) = aip(x, xo,s), 5ip(x, xo,s) = — aip(x,xo,s). (3.53) 

This is a symmetry of the folded theory discussed in appendix B. It is obviously inap- 
propriate to break this vector U(l) symmetry in the bulk. The "chiral" symmetry at the 
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boundary can be broken by a physical quark mass term mfqq keeping the bulk vector 
symmetry. But this is forbidden by the chiral Ward-Takahashi identity in the bulk. The 
symmetry should be broken only at the boundary. 

We introduce boundary source fields as a component of projected out degrees of free- 
dom in (3.18) and (3.19) 

X(x,s) = P^{x,0,s), X'(x,s) = P + ip(x,N T ,s), (3.54) 

X(x, s) = (^ P_) (x, 0, s), X'(x, s) = (^ P+) (x, N T , s) (3.55) 

The orbifolding symmetry breaking term takes the form 

^breaking = X(x, s)O st P+1p(x,0,t) + (ip P+) (x, 0, s)O st X(x, t) 

+x'(x,s)d s tP^(x,N T ,t) + (^P_) {x,N T ,s)6 st X'{x,t), (3.56) 

where O is a local operator which anti-commute with T = 7o75-PQ- Candidates of O are 
70, 75, Q, P and 

K(u) st = (P l 5 s ,n 5 +i-u + Pr$s,u) {Pl5i,u + PRSt,N 5 +i-u) , (3.57) 

K(u) st = (PiA,JV 5 +i-u + P L S s ,u) (Pn$t,u + PLSt,N 6 +i-u) , (3-58) 

where K(l) produces a physical quark mass term t/jK(l)ip = qq. 

Among these candidates 75 and Q are forbidden by the parity symmetry. 70 is not 
consistent with the charge conjugation. P, K{u) and K(u) break chiral symmetry, which 
however is broken at the boundary. Since P, K(u) and K(u) are consistent with parity 
and charge conjugation they are proper candidates of orbifolding symmetry breaking term. 
However we adopt only K(l) in this paper to reproduce the same surface term as the 
continuum theory. We define the surface term as 



Surface = -a 3 ^(\\(x,s)^(l) st P+V>0?,O,£) + (i/, P+) (x,0, s)K(l) st X(x,t) 

X 

+x'(x,s)K(l) st P^(x,N T ,t) + (^P_) (x,N T ,s)K(l) st X\x,ty 

J2{- p(x)P-q{x)\ X0=0 ~ q( x ) p +p(x)\ xo =o 

X 

p'(x)P + q(x)\ xo=NT - q(x)P^p\x)\ xo=NT ), (3.59) 



a 3 ' 



where q and q are active dynamical fields at the temporal boundary, p and p are boundary 
source fields for the physical quark fields 



P+q(x)\ xo=0 = p(x), P-q(x)\ X0=NT = p'(x), (3.60) 

q(x)P-\ X0=0 = p(x), q(x)P+\ X0 =N T =p'(x). (3.61) 



Since this surface term is not a general orbifolding symmetry breaking term, general 
anticipation is that we will need all sort of breaking terms to renormalize quantum correc- 
tions. However if we consider Green functions constructed with physical quark operators 
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only we may expect that quantum corrections which appear in these Green functions can 
be renormalized into a shift of physical operators and physical quark source fields p, ~p, 
p' and ~p' . The situation is similar to the physical quark mass term mfqq. This term is 
not a general chiral symmetry breaking term. However if we consider a Green function 
of the physical quark fields only, quantum corrections which appear with chiral symmetry 
breaking can be renormalized into the physical operators and mass term. We did not need 
all the breaking terms for renormalization if we consider the physical quark Green function. 
Explicit calculation is necessary to confirm this expectation for source fields. 

We check validity of this surface term at tree level. According to Ref. [10] we introduce 
the generating functional 



z f [p',p';p,p;v,V,U] = DipDi)e^p^-S F [U,ip,^;p',p',p,p] 

+a 4 2J {ip(x, s)rj(x, s) + r](x, s)ip(x, s)) > , (3.62) 

x,s ' 

where r){x) and rj(x) are source fields for the fermion fields and the total action Sf is given 
as a sum of the bulk action (3.26) and the surface term (3.59). We notice that the fermion 
fields ip and ifj obey the orbifolding condition (3.17). The correlation functions between 
the boundary fields are derived with the same procedure as Ref. [10]. 

(il>(x, s)^(y, t)) = GZi(x, y; s, t), (3.63) 

(q(x)q(y)) = Gl* axk (x,y), (3.64) 

(q(x)((y)) = Gl^ Tk (x,y)P + \ yo=Q (3.65) 

(<z(*)C'(y)) = GqLk(*, y)P- \ yo=NT , (3-66) 

(C(x)q(y)) = P- Gll ark (x,y)\ xo=0 , (3.67) 

(C'(x)q(y)) = P+ G^ aik (x,y)\ xo=NT , (3.68) 

{((my)) = P- Gll ark (x,y)P+\ X0=0m=0 , (3.69) 

(C(x)C'(y)) = P GH ark (x, y)P^ \ Xo=0:yo=Nr , (3.70) 

(C'(x)((y)) = P + Gll lk (x,y)P + \ XQ=NTM=0 , (3.71) 

(C(x)c'(y)) = p + Gll lk (x, y)P-\ X0=NTM=NT ■ (3.72) 

The propagator G^ f and G^ ark are given in (3.31) and (3.38). We notice that the above 
propagators between the boundary fields and physical quark fields approach to the contin- 
uum SF boundary propagator without any 0(a) term at tree level. 

3.4 Effective action of the domain-wall fermion 

In order to perform numerical simulation with dynamical fermion we need to introduce the 
Pauli-Villars field to cancel bulk contribution in fifth direction. The Pauli-Villars field is a 
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four component complex scalar and its action is given by 

N T N 5 

Spy = ^2 ^2 5Z ~$( x > s ) D Pv(x,v,s,t)<i>(y,t), ( 3 - 73 ) 

x,y xo,yo=—NT+i s,t=i 

where Dirac operator for the Pauli-Villars field is given in the same form as the domain-wall 
fermion Dirac operator (2.6) with rrif = 1 

^pv = D dwi (m f = 1). (3.74) 

This Dirac operator does not commute with the orbifolding operator A = 7075 PQR because 
of the mass term. It is not straightforward to introduce the Pauli-Villars field by orbifolding. 
In this paper we propose to implement it by the effective Dirac operator [34, 35]. 

The effective Dirac operator appears in an effective action of the physical quark field 
(2.9) (2.10) and "physical" Pauli-Villars field 

Q{x) = (P L S a ,i + Pr8,,n 5 ) <!>{x, s), (3.75) 

Q(x) = 4(x, s) (5 S , N5 P L + S Stl P R ) . (3.76) 

The effective action is given by integrating out all the bulk fields other than physical fields 
at the fifth dimensional boundary [34] 

Seff = q(x) (D cS ) xy q(y) + Q(x) (D cS + l) xy Q(y). (3.77) 

In its derivation the effective Dirac operator D e g is given as an inverse of the full physical 
quark propagator 

D cS = -^-, (3.78) 



(qq) 



whose explicit form is 



1 + 75,5 i-T^ 5 1-H' , 1 , 

Here Dw is a four dimensional Wilson Dirac operator with negative mass —2 < —M < 0. 
In N$ — ► 00 limit the Dirac operator becomes 



D cS = >^f, T = e~ H , (3.80) 

l-7se [H 



where e(x) is a sign function 

e(x) = -£=. (3.81) 

We can easily check that this Dirac operator is exactly chiral symmetric 

{75,^cff} = (3.82) 
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in N§ — ► oo limit and should be non-local to satisfy the Nielsen-Ninomiya's no-go theorem. 
The effective Dirac operator is related to the original domain-wall fermion and the 
Pauli-Villars Dirac operator through determinant 

det — — D dwi = det - c = det D Ns , (3.83) 

-L>pv JJ eS + 1 

where Dn 5 is a truncated overlap Dirac operator [34, 35]. Hereafter we take N$ — ► 00 limit 
implicitly and write Z?_/v 5 ^oo = -Cod- In terms of the domain- wall fermion the overlap 
Dirac operator is defined as 

D 0» = j^~. (3.84) 

£>OD satisfies the Ginsparg- Wilson relation [36] 

{75, Dod} = 2L>od75 A)D- (3.85) 

If we introduce physical quark mass term we have a massive overlap Dirac operator through 
determinant 

A?D(m/)= n ff + ? 7 = -Pod +m f (l-D D) ■ (3.86) 

The effective Dirac operator of the orbifolded domain-wall fermion system is defined 
in a similar way. Since the four dimensional Wilson Dirac operator Dw commute with the 
four dimensional time reflection operator S = i^^qR we have following anti-commutation 
relations 

{£,#'}=(), |S,#}=0. (3.87) 

By using these relations we can easily show that the effective Dirac operator (3.80) anti- 
commute with the four dimensional orbifolding operator Y defined in (3.23) 

{Y,D cS } = 0. (3.88) 

The massless overlap Dirac operator defined in the above satisfy "Ginsparg- Wilson rela- 
tion" for the orbifolding transformation [25] 

{r,D D} = 2r>0Drr»0D. (3.89) 

We define the Schrodinger functional effective Dirac operator as an inverse of the 
orbifolded full quark propagator (3.38) 

D e S | = H+ -4yIL = U + D cS U , (3.90) 

where inverse means that in a sub-space 

^cflfCquark = 211+) G quark D cfl = 211^. (3.91) 
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Contribution from the Pauli-Villars field is introduced to reproduce the Schrodinger func- 
tional overlap Dirac operator defined in Ref. [25] § 

z>gb = n +^eff n_ - 1 = n+Dooii- , (3.92) 



where 



n ± = i^I, r = r(l-2£>oD). (3-93) 



This is not a unique definition of the SF overlap Dirac operator but we can define 
another Dirac operator as 

^od = jr^—U+D^U.. (3.94) 

JJ e s + i 

These two Dirac operators are related by unitary operators 

« = — ^— (1 - 2Z?od) + —£-, u' = 75^75 (3.95) 

by 

uD^ = DZ, ^D^u' = D S J d . (3.96) 

Here we used a fact that the effective and the overlap Dirac operators commute with the 
four dimensional time reflection operator S. 

As was discussed in Ref. [25] the SF overlap Dirac operator does not have 75 Hermiticity 
relation. Instead we have 

(Z?S F D ) f = 75^0075- (3.97) 

In order to define real fermion determinant we may need even numbers of flavours and 
different Dirac operators for each flavours. An example for two flavours case is 

/ nSF \ 

4? = ( ° D nSF)- (3-98) 

We notice that U(2) vector flavour symmetry is broken to U(l) x U(l). Determinant of 
this Dirac operator is 

det 4? = det 4? 75 = det ( U + D cS 1 * d\ ff II + J , (3.99) 

which is re-written in terms of pseudo-fermion field x 

,(2) 



det D K SF 



jv (n +x f ) v(u +x )cx P f- x f n + (-±- + 1 j (-±- + 1) u +x ) (3.100) 



^The Ginsparg- Wilson relation (3.89) in this paper is different from that in Ref. [25] by factor two and 
so is the definition of F. 
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The determinant is denned in a sub-space 1i- = {^|n_?/> = 0} of eigenfunctions. In eval- 
uation of the fermion force we need to calculate 

-^ + l) = «W> + I)" 1 , (3-101) 

which corresponds to inverse of the overlap Dirac operator. 

The orbifolded effective Dirac operator is modified as follows when we introduce the 
mass term (3.29) 

DH(m f ) = ^II + (D cS + m fV ) II_. (3.102) 

Taking into account a contribution from the Pauli-Villars Dirac operator the massive SF 
overlap Dirac operator is defined as 

£>od(™/) = \^+ (Afl + m f n) n -^7T7 = \^+ Pod + m f n (1 - D OD )) U-, 

(3.103) 

op 1 1 

Dov(m f ) = 2^^T n + (Afl + m fV) n - ( 3 - 104 ) 

Although we do not have a unitary transformation to relate D^^^mf) and -D D( m /) we 
have a Hermiticity relation 

(^OdC™/))* = 75^od(™/)75. (3.105) 

We also need even numbers of flavours to define a real fermion determinant. 

4. Schrodinger functional with twisted boundary condition 

In the previous section we presented an orbifolding formulation of domain-wall fermion in 
finite box, in which the homogeneous proper boundary condition (3.24) (3.25) is satisfied. 
This is a solution of our purpose to define a finite volume renormalization scheme. However 
this may not be the unique solution of our requirement that the theory has a mass gap and 
is kept to be renormalizable in a finite box. In this section we propose another orbifolding 
formulation to adopt chirally twisted boundary condition [25, 37, 38]. As was discussed in 
Ref. [25] the chirally twisted boundary condition has advantages that the fermion deter- 
minant becomes real and the mass term is introduced easier. For domain- wall fermion the 
Pauli-Villars field can be treated in a straightforward way by orbifolding. 

4.1 Orbifolding construction of chirally twisted boundary condition 

As will be discussed later the fermion determinant becomes real for even numbers of 
flavours. In this section we adopt two flavours case for instance. We start from the 
massless orbifolded action (3.26) and introduce the twisted orbifolding by chirally rotating 
the fermion field 

^e^f^y, lp = ^'e-^ QT \ (4.1) 
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where r 3 is the Pauli matrix to act on flavour space and Q is the vector charge (3.8) for 
chiral transformation. In terms of the rotated field the orbifolded action is given by 

Ssf = J2 ^' 5 dwfV>', DZf = n_D dwf n_, (4.2) 

where 

is a twisted orbifolding projection with time reflection operator £ defined in (3.1). 

As was discussed in sub-section 3.1 the Dirac operator has no index and the chiral 
transformation is not anomalous even for Abelian case. This formulation with twisted 
orbifolding projection is equivalent to the original one for massless theory. Here we notice 
that the twisted orbifolding operator Sr 3 commute with the massive domain-wall fermion 
Dirac operator 

[Sr 3 ,D dwf (m/)] =0 (4.4) 

since we adopted time reflection invariant gauge configuration. We can extend this twisted 
formulation to massive theory 

S'Zf = E ^ 5 "wf (m/M ^dwf = 5-Aiwf (m/)5_. (4.5) 

It is straightforward to introduce the Pauli- Villars field through orbifolding 

«5pv st = E^pv^ 5f£ = n_i>pvn_ (4.6) 

since .Dpy = -Ddwf(m/ = 1) and is commutable with the orbifolding operator. 
The fermion fields satisfy the twisted orbifolding projection condition 

n+v = o, ?n + = o, (4.7) 

which brings the following boundary conditions 

P+if)(x,Q,8) = 0, P-ip(x,N T ,s) = (4.8) 

(^P + )(x,0,a) = 0, (^pJ)(x,N t ,s)=0 (4.9) 

with projection operator 

P ± = l± ^oPr\ (4 1Q) 

In terms of the physical quark field the projection condition becomes 

(P L 5 S>1 + P R 5 s , Ns ) (U + ) ij(x,t) = U+q(x) = 0, (4.11) 

V / s,t 

$(x, t) (n + ) (S s ,n 5 Pl + 8s,iPr) = q(x)U+ = 0, (4.12) 

V / t,s 

n± = X± Sr , S = i 7 57oP, (4.13) 
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where £ is the time reflection operator in four dimensions. The boundary condition for 
the physical quark field is 

P + q{x)\ X0=0 = 0, P-q(x)\ X0=NT = 0, (4.14) 

q(x)P+\ X0=0 = 0, q(x)P-\ X0=NT =0, (4.15) 

P ± = 1± ^ 57or3 . (4.16) 

We have two comments. The orbifolded Dirac operator with twisted projection has a 
following Hermiticity relation 

D SF (m? = ^T^DsFim^r 1 ' 2 , (4.17) 

which is also the same for the orbifolded Pauli-Villars Dirac operator. The SU{2) flavour 
symmetry is broken to U{\)y x U{\)?, as in the chirally twisted mass QCD. 

4.2 Free propagator 

The original theory before orbifolding has a mass gap proportional to 1/T because of anti- 
periodicity in temporal direction. This property is robust against orbifolding process and 
survive in the twisted orbifolding formulation. We will check this property at tree level by 
using propagator. 

The fermion propagator is defined as an inverse of the orbifolded Dirac operator in a 
sub-space 

GZ { (x,y,s,t) = 2(D s d l f Y l = 2 (rL-i-rL^ , (4.is) 



x,y;s,t 



i — z, i ii_ 

/ x,y;s,t \ -Ddwf 

«f [DZl) = {DZ^' 1 DZ = S_. (4.19) 

At tree level this propagator can be written in a simple form as 

G<gt(x, y;M) = ^E e^^GZM so, W, s, t), (4.20) 



N? 



G¥Av-xn un-s t) = - V" I - I i ( e ip ^ Xo ~ Vo "> - e ip °( Xo+yo A (Pa.~) 

G dwf{P ,x ,yo, S ,t) 2aNT ^2.^ \D M (p))^\\ e ) V + )t>,t 

( e ipo(xo-ya) _|_ e ipo(xa+ya)\ (p_\ I. (4.21 



+ V 

-Cdwf (p) is the domain- wall fermion Dirac operator in momentum space without orbifolding 
projection, whose inverse is given in appendix C. We notice that the temporal momentum 
Pq satisfies the quantization condition (3.35) and there is no extra fermion zero mode. 

The physical quark propagator is given by selecting the contribution from the boundary 
fields in fifth direction 

<^Lk(s, y) = (PlS s ,i + PrS s ,n 5 ) GZt (s, y; a, t) (6 t , N5 P L + 5 ttl P R ) 

= 2 (fl_G quark n_) , (4.22) 

v / x,y 
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where G quar k(x,y) is the physical quark propagator in 2Nt x N^ space-time without any 
projection. Following Dirichlet boundary conditions 

P + d® uttk (x,y)\ X0=o = 0, P-G^ k (x,y)\ X0=NT = 0, (4.23) 

G^Lkfo y)\yo=oh = 0, G^ ark ( x , y)\ yo=NT P- = (4.24) 

are satisfied for this quark propagator. By ignoring sub-leading terms in e 5 the propa- 
gator takes the following form at tree level 



Z_j quarkV > U) ^^ ^ \l - e a W (p) J 



n=-N T +l v Ky ' ■ 

x | ( e - ipo ^ - e ipoyo ) P+ + ( e - ipoyo + e ipoyo ) P_} (4.25) 

We emphasize that the physical quark has a gap proportional to 1/T because of the anti- 
periodicity (3.35). This formulation satisfies one of the requirement. 

4.3 Surface term 

In this subsection we consider a twisted orbifolding symmetry and introduce a coupling 
to the boundary source field (surface term) as a symmetry breaking term. The orbifolded 
action (4.5) is invariant under the following twisted orbifolding transformation 

5(U-i?j(x,s) = a(U-i?j(x,s), <$te_)(x,s) = -aft/fi_)(a;,s), (4.26) 

where remaining degrees of freedom TL+ip and Y>n + are kept intact. The boundary source 
fields are elements of H+ip and ipH+. 

We define a surface term as an orbifolding symmetry breaking term, which is consistent 
with parity 

ip(x, s) -»■ 7oP s tT 1 'V(-x, x ,t), tp(x, s) -»■ ip(-x, xq, i)7 P ts r 1 ' 2 , (4.27) 

charge conjugation 

if>(x, s) - CP st r 1 ' 2 ^ T (x, £), ^(x, a) - t/; t (x, t) (-C" 1 ) P^r 1 ' 2 , C = 727 o (4.28) 

and 

tfj(x,s)^C(PQ) st ^ T (x,t), ^(x,s)^ip T (x,t)(-C~ 1 )(PQ) ts , C = 7 27o (4.29) 

and vector U(l)3 symmetry 

Sip(x, s) = (3r 3 ij(x, s), 5lp(x, s) = -,#(x, a)r 3 . (4.30) 

of the orbifolded domain-wall fermion. Here we modified the parity and the charge conju- 
gation transformation to be consistent with the twisted orbifolding projection. 
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Using the orbifolding projection (4.7) the orbifolding transformation (4.26) is shown 
to be a "chiral" transformation at the boundary in which a half of degrees is rotated 

5 (P-A (x, 0,s)=a (P-A (x, 0, s), (4.31) 

5 fe_) (x, 0, s) = -a fep_) (x, 0, s), (4.32) 

5 (P+V) (x, N T , s) = a (p+^j (x, N T , s), (4.33) 

8 (^P+) (x, N T ,s) = -a (^P + ) (x, N T , s) (4.34) 
and is a vector U(l) transformation in the bulk < xq < Nt 

5ip(x, xo, s) = aip(x, xq, s), 5ip(x, xo, s) = —aip(x, xo, s). (4.35) 

The symmetry should be broken at the boundary. 

We introduce boundary source fields as a component of projected out degrees of free- 
dom 

X(x,s) = P + ip(x,0,s), \'(x,s) = P-ip(x,N T ,s), (4.36) 

X(x, s) = (^ P+) (x , 0, s), \'(x, s) = ($ P_) (x, N T , s). (4.37) 

The orbifolding symmetry breaking term takes the form 

^breaking = \(x, s)d st P-1p(x,0,t) + \\j) P-J (x , , s)0 st X(x , t) 

+\\x,s)d st P + 4>(x,N T ,t) + (^P+) (x,N T ,s)d st \'{x,t), (4.38) 

where O is a local operator which anti-commute with i^^qPt 3 . Candidates of O which is 
consistent with the parity, charge conjugation and ^(1)3 symmetries are PQt 3 , K(u)Qt 3 
and K(u)Qt 3 , where K and K are defined in (3.57) (3.58). 

As in the previous section we adopt K(1)Q for the surface term to couple only to the 
physical quark field 

Surface = ~a 3 ^ fe s)(K(l)Q) st T 3 P_ ^(x , 0, t) + fop_) (x, 0, s)(K(l)Q) st T 3 \(x , t) 



x =0 



+A {x, s){K(l)Q) st T 3 P + ^{x, N T , t) + (V P+) (x, N T , s)(K(l)Q) st r 3 X'(x, t) ' 

= a 3 V (- p(x)j 5 T 3 P_q(x) - q(x)P^ 5 T 3 p(x) 

^"T ^ 2:0=0 

x 

- 7f(x) 75 r 3 P + q(x) -q(x)P +lb T 3 p'(x) ). (4.39) 

p and ~p are boundary source fields for the physical quark fields 

P+q(x) U =o = P(x), P-q(x)\x =N T = p'{x), (4.40) 

q(x)P+\ X0=0 = p(x), q{x)P-\x =N T =p'(x). (4.41) 

Although this surface term is not a general symmetry breaking term, we also expect 
that quantum corrections can be renormalized into a shift of physical operators and physical 
quark source fields p, ~p, p' and ~p' if we consider Green functions constructed with physical 
quark operators only. 
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4.4 Effective action of the domain-wall fermion 



For the twisted orbifolding formulation of finite volume field theory the Pauli-Villars field 
is introduced directly as in (4.6). Total contributions from fermion and Pauli-Villars field 

is 



detn 



Dpv 



-Ddwf(m/)n_ 



(4.42) 



where the determinant is defined in a sub-space TI+ = < iplU+t/j = > of eigenfunctions. 
In this sub-section we will show that this determinant is equivalent to that of the overlap 
Dirac operator with twisted orbifolding [25] 

detfL— — £> dwf (m/)fL = detfL£> D(m/)n_. (4.43) 



n+ 



Dpv 



m 



For this purpose we adopt the Schur decomposition procedure for the effective Dirac 
operator [39, 40]. Statement of the Schur decomposition is that the overlap Dirac operator 
is given as a Schur complement of the domain-wall fermion Dirac operator divided by the 
Pauli-Villars Dirac operator 



1 



D PV 



-D, 



'dwf(«VJ 



VU- l (l)D { ^(m f )U(m f )Vl 



(4.44) 



,(5) 



Here V, U(rrif) and -Dqq(?71j) are matrices in fifth dimension and their explicit forms for 
N$ = 6 case are given by 



V 



U(m f ) 



(Pr 






Pl\ 




Pl Pr 








Pl Pr 






= p R + n 


Pl 


Pr 








Pl 


Pr 




\ 




Pl Pr) 






(I 




-T (P L - m/Pfl) \ 




1 




-T 2 (P L - m f P R ) 


\ — 


1 




-T 3 {P L - m f P R ) 


) — 




1 


-T 4 (P L - m f P R ) 






1 


-T 5 {P L 


- mfP R ) 



■1)Pl, 



(4.45) 



(4.46) 



/I 



D ( ol(m f ) 



\ 



(4.47) 



where .DoD(m/) is a truncated four dimensional massive overlap Dirac operator 
A>d(to/) = \ (1 + 75^) + m f ( 1 - ^ (1 + 75^) 



(4.48) 
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with the same definition for S in (3.79). The truncated Dirac operator turns out to be the 
ordinary overlap Dirac operator (3.86) in TVs — ► oo limit. Q~(rrif) is a hopping operator in 
fifth direction (2.8). So we have 



det — — D dw f(mt) = det D OD (mf) 

JJpy 



(4.49) 



for ordinary domain-wall fermion system. 

We start from the orbifolded domain-wall fermion Dirac operator divided by the Pauli- 
Villars Dirac operator 



~ 1 



D^ = n_-— D dwf (m/)n_ = U_VU-\l)D^(m f )U(m f )P^U_. 
Upv 



(4.50) 



We consider multiplication of the projection operator on unitary matrix V and we have 

(4.51) 



~ „, ~ m £ £ l±Pfi-(-l)Er 3 

n_p = n_pn_, n± = - — - — . 



We notice that a matrix in the projection II ± has a following form 

where P(n 6 -i) 1S a (-^5 — 1) x (-^5 ~~ 1) matrix of the form 

/0 1\ 
10 



Pi 



(iVs-l) 



10 

10 

\\ 0/ 



(iV 5 = 6). 



The projection operator II-|- is written as a direct sum of two projections 



n± 



n± 



n± 



where 



ll-i- — ; 



(4.52) 



(4.53) 



(4.54) 



(4.55) 



is a projection operator in N$ — 1 sub-space. 

Taking into account the explicit form of the matrix U(m.f) its determinant multiplied 
by the projection becomes 



det U{m)Ii. 

(+subspace) 



det U-U-Hm) 



(+subspacc) 



det n_C/(m)n_ = det \ U „ j = 1, (4.56) 

(+subspace) (+subspace) \ I I _ / 

_ __ ^ /~(iV 5 -l) \ 

= det n_t/~ 1 (m)n_ = det n - _ ) = 1. 

(+subspace) (+subspace) \ 1 I _ / 

(4.57) 
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Substituting this relation determinant of the total Dirac operator is equivalent to that of 
the orbifolded overlap Dirac operator 

det-D^F = detft_pn_C/- 1 (l)n_ J Dgi(m / )n_C/(m / )n_'P+n_ 

= ^ dfit (ii)8( m/ )i 

(+subspacc) \ 

/-(JVs-l) \ /1 N /~(JV 5 -1) \ 

= det f n - ) f 1 ^" 1 ) ) U ) 

(+subs P acc) y n_y V ^odK)/ ^ n_/ 

= detn_DoD("i/)fl„ (4.58) 

and we get expected result. 

At last we have a comment on Hermiticity. The five dimensional total Dirac operator 

(5) 

Z)g F has a following Hermiticity relation 

D$ = Tst^fV 1 ' 2 (4-59) 

and its determinant is real. Since our domain-wall fermion Dirac operator does not have 
index the chiral rotation (4.1) is well defined even for single flavour case and we can define 
a single flavour orbifolded Dirac operator as 

^SF SlC = ^s^Wm/)^. (4.60) 

However we do not have a Hermiticity relation for this Dirac operator and the determinant 
is not shown to be real. We may need even numbers of flavours to avoid this problem. 



5. Conclusion 

In this paper the orbifolding formulation of the finite volume field theory is applied to the 
domain-wall fermion. In order to reproduce the proper SF Dirichlet boundary condition 
we need both the time reflection and the chiral symmetries. Application of this procedure 
to fermionic part is straightforward because of good chiral symmetry of the domain-wall 
fermion. Since there is no chiral symmetry for the Pauli-Villars field it is introduced by 
using the effective Dirac operator to reproduce the SF overlap Dirac operator. The surface 
term is given as an external source field to break the orbifolding symmetry. 

The SF Dirichlet boundary condition may not be the unique choice to define a finite 
volume field theory suitable for renormalization scheme. A finite volume field theory with 
chirally twisted boundary condition is also proposed. Time reflection symmetry is enough 
to reproduce the twisted boundary condition by orbifolding. We can treat the fermionic 
part and the Pauli-Villars field in an equal footing. We have a 75 Hermiticity relation 
for the orbifolded Dirac operator and the total determinant is real. This formulation is 
applicable to two flavours dynamical simulation. 
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A. Effective action of chiral symmetric Dirac operator 

In this appendix we derive the effective Dirac operator of the physical quark field for an 
action with the chiral symmetric Dirac operator (3.12). Four dimensional part of the 
symmetric Dirac operator is the same as the ordinary Dirac operator (2.6). Hopping term 
of this Dirac operator into the fifth direction takes the form 



P L n + (m f = 0) + PflJT (m/ = 0) 



( ° 


Pl 











\ 


Pr 





Pl 














Pr 


























Pl 














Pr 





Pl 


V o 











Pr 


/ 



(A.l) 



for massless case. If there were no quark mass this Dirac operator is equivalent to two 
independent domain-wall fermion with half fifth dimensional size of N$/2. It is easily 
shown that there are two extra zero mode (doublers) at the middle boundary s = -^- and 
s = -£■ + 1 related to the exact chiral symmetry at finite N$. 

The physical quark fields may be defined in the same manner as (2.9) and (2.10). We 
can integrate out the bulk field other than q and q according to Ref. [34, 35, 39]. We start 
by writing the fermion field as a vector in fifth direction and chirality. For N§ = 6 we have 



tf J 



(-01H i>lL 4>2R 1p2L 1p3R 4>3L ^AR V>4L 4>5R V>5L V>6R V>6L ) 



* = ( tPlL ll>lR i>2L ^2R ^3L ^3R 4 AL 4 'AR Y^L ^5R 4&L 4&R ) , 



where 



i>R/L = Pr/l4, 4r/l = ^Pl/r- 



Then we change variable as 



* 



>T 



(iplL IplR i>2L 1p2R i>3L ^3R ^AL 4>AR ^5L 4>hR V>6L 



D, 



( a (3 

a (3 
0o «o 



dwf 



where 



a 



B -C* 

-1 



«o = Prcx, 



\ 



a (3 

a (3 
fa a J 



(A.2) 



* = ( V'lH 4>2L i>2R i>3L ^3R $1L i 'W ^5L V>5R 4&L Ar 
The Dirac operator is written as follows in terms of the primed field 



foil), 

4al)- 



(A.3) 



(A.4) 



25 



1 ( c B ) h -1'iJ- (A.o) 



C. 
B 



Xjl 



°>2 (<W,s/£ / m( x ) ~ 5 x-n, y ul{yyj , (A.6) 

(1 - M)6 xy - - (Sx+wUnix) + 5 x .^ y Ul(y) - 2S xy ^ , (A.7) 



(7 , 



7,= ( t H- (Ai 



We integrate out all the fields except for the physical quark field 

q{x) = P L 4>(x, 1) + P R il>{x, 7V 5 ) = Pl7o^(x, 1) + Prj ^'(x, N 5 ), (A.9) 

q(x) = ^(x, l)P R + lp(x, N 5 )P L = ^'(x, ^-)j Pr + ?(*, /V 5 )7o-Pl (A.10) 

according to Ref . [34] . Result is given as a full quark propagator 

^ )= Ki~w 5 ) = ^- (A - n) 

Here D e g is the truncated effective Dirac operator (3.79) with half size of fifth dimensional 
length 

1 - 75^ i + y 2 

Transfer matrix is given by 

T = 7570 (-a^ 1 ) 7075 = ^-^7- (A.13) 

The full quark propagator anti-commutes with 75 even at finite N$. In A5 — > 00 limit the 
effective Dirac operator D e g anti-commutes with 75 exactly and the effective Dirac operator 
-^eff™ w ith symmetric construction becomes the same as that of the ordinary domain-wall 
fermion D e g. 

We introduce the Pauli-Villars field in the same manner with the Dirac operator 

£>p y v = Aiwf (m/ = 1) + X. (A.14) 

The effective action of the physical quark field q, q and the physical Pauli-Villars field Q, 
Q is given by 

S cS = qD% m q + Q {D s ^ m + l) Q. (A.15) 

The overlap Dirac operator is given to reproduce the same determinant as the effective 
action 

D OB = n sym T -, ■ (A. 16) 
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Because of exact chiral symmetry of D^ 1 the overlap Dirac operator Dq^ 1 satisfies the 
Ginsparg- Wilson relation even at finite N§. 

Compensation of the exact chirality at finite ./V5 is a non-locality in the overlap Dirac 
operator, which comes from the extra zero mode in the middle of fifth direction. However 
we can show that the non-locality is exponentially small in N§ and disappears in N5 — ► 00 
limit. In order to extract the non-locality we define explicit breaking term of the chiral 
symmetry of the ordinary effective Dirac operator (3.79) or the truncated overlap Dirac 
operator at finite N5 

6n > = 75 i^k) + i^k) 75 = 15 D^m + ^dW 5 - 275 - (A - 17) 

The chiral symmetric effective Dirac operator is re-written as 

1 II 

7«m= 7, T^iVs, ( A - 18 ) 

DJs D cS 2 2 

where we used a fact that the breaking term commutes with 75 

[S N5 , l5 ]=0. (A.19) 

The chiral symmetric overlap Dirac operator is given in a following form 

^o y D = : ,' -. %■ (A.20) 

1 - 2-^0075^^ 

Dob m denominator may bring a non-local factor into the overlap Dirac operator. However 
as was shown in Ref. [34] 5n 5 is exponentially small in N§. The physical part of the chiral 
symmetric Dirac operator -D^wf c °i nc ides with that of the ordinary Dirac operator in 
N$ — ► 00 limit. 

B. Folding of temporal direction 

In our formulation with the orbifolding (3.17) fermion fields in negative time — Nf < xq < 
can be written in term of those in the positive region 

ip(x,-x ,s) = (T) st tp(x,x ,t), r = 7o7 5j PQ- (B.l) 

Half of the field degrees of freedom can be eliminated explicitly by folding the temporal 
axis into the non-negative range < xq < Nt together with the boundary condition (3.18) 
(3.19). 

For this purpose we introduce four projection operators in temporal direction 

T_ for -N T + 1 < x < -1, 

To for Xq = 0, 

T+ for 1 < x < N T - 1, 

Tt for xq = Nt, 
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which pick up the fermion fields on the corresponding region. For instance 



i-r , , a J V'(^o) for 1 < x < N T - 1 

(/+) — ^°) = { otherwise " (&2) 



'xo,yo 



Summing up four projection operators we have a unity 

1 = T_ + T + T+ + T T (B.3) 

and T a 's have a projection property 

T a T p = T a 5 a:f3 . (B.4) 

These projection operators satisfy the following relation with the time reflection operator 
R 

RT + = T-R, i?T_ = T + R, RT = T R = T , RT T = T T R = -T T . (B.5) 

By using these properties we have an identity relation 

n+ = n+ (T+ + T_ + T + T T ) 

= n+ (T+ + T + T T ) (2T+ + T + T T ) U + (B.6) 

and the orbifolded action (3.26) can be re-written in terms of the fermion fields depending 
on the non-negative region only 

N T 

S = I2 E ^'(x,s)D^ cd (x, y] s,t)^\y,t), ( B - 7 ) 

x,y xo,yo=0 s,t 

where ip" and ip are defined as 

Y>"(x, xo, s) = ((T + + T + T T ) TL + ) xoyo . st 4>(x, yo, t), (B.8) 

^"{x, x , s) = lf(x, y , t) (n+ (T+ + T + T T )) yo ^. ts , (B.9) 

which have no dependence on negative time. These fields can further be written as 



r(x,x ,s) = (T + + T P + + TTP-) Xoyo . s ^(x,yo,t), (B.10) 

?'(x,x 0) s) =^(x,y ,t)(T + + P + T + P_T T ) yoxo . As (B.ll) 



by using (B.5) and identification (B.l). There is no constraint on positive bulk fields. 

jSF 



The folded Dirac operator D^L f is given formally as 



D^f d = \ (2T + + T + T T ) n + D dwf n + (2T+ + T + T T ) . (B.12) 
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This Dirac operator can be written in more explicit form by using the orbifolding symmetry 
(3.16) and the ultra local property of the domain-wall fermion Dirac operator, with which 
we eliminate the term like T + D dw fAT + = T + D dw fT-A 

Aiwf C = 7;ToP+DdwfP+To + T P + Dd w {T + + T+.Ddwf-P+^b + T + D dw {T + 

+T T P^D dw{ T + + T+D dw{ P^T T + ^T T P^D dwi P„T T (B.13) 

(p + !Bp_p + -P +P _ Uo{0) \ 

-P+P + E/J(0) D( 3+1 ) -P-Uo(l) 

-P+U^l) Z)( 3+1 ) -P-U (2) 

-P + Ul{2) Z)( 3+1 ) -P_P_ C/ (3) 

V -P-pM$) P-^P-P-J 

(B.14) 

where the matrix represents the Dirac operator in temporal direction for Nt = 4. Z)( 3+1 ) 
is the Dirac operator in spatial direction and the fifth direction 

D {3+1) {x,y;s,t) = ( — % -U i (x)S yiiXi+1 + — -UJ {y)8 yuXi - 1 J S XOtVo 8 a> t 

+ {^f^^tt + Z ^P^) S *,v + ( 5 " M)8 x , y 8 a , t . (B.15) 

There is no constraint for the bulk region 1 < Xo,yo < Nt — 1, which is nothing but 
ordinary domain-wall fermion Dirac operator. 

We notice that the projection operator P± at the boundary does not commute with 
the 70 chiral projection P±. If we consider an eigenvalue equation of this Dirac operator a 
zero mode dumping solution 

ip = P_(l - M) x ° + P+(l - M)( Nt - x °) (B.16) 

in temporal direction, which have broken the chiral symmetry "dynamically" in a naive 
formulation, is forbidden by this boundary term. 

The fermion propagator is given as an inverse of the folded Dirac operator 

G™ = 2 (T+ + T + T T ) n +J D-^ f n + (T + + T + T T ) , (B.17) 

where the inverse is defined in the ordinary meaning for the positive bulk region < xq < 
Nt and in terms of the projected sub-space at the boundary 

^dwf ed G^ d f ed = T+ + P+T + P^T T . (B.18) 

C. Free fermion propagator 

Inverse of the massless domain- wall fermion Dirac operator in momentum space is derived 
according to the procedure of Ref. [30]. In this appendix we omit derivation and give the 
result: 

1 



T>dwf(p) 



-i^sinp^ + W-Q-) G R P L + (-^sin^ + W - Q+) G L P R , (C.l) 
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where £1 and W are denned in (2.8) and (3.37). G R and Gl are denned as 

G R (s, t) = G°(s-t) + A ++ e a{s+t ^ + A+^e^-V + A_ + e a( -- s+ ^ + i__e a (" s ^(C.2) 
G L (s,t) = G°(s - t) + B ++ e a ^ s+t ^ + B + ^e a{s ~^ + B^ + e a( - s+t ^ + B^e^" 3 ^ (C.3) 
G°(s -t) = C ( e <*W-\'-t\) + e -a(JV 8 -|5-t|)^ ( C _ 4 ) 

with exponent and coefficients given by 

l + W 2 +sm 2 Pfl 
cosha = — £, (C.5) 

(^ (C.6) 

4VFsinh a sinh (aiV^)' 

A ++ = F(l-We- a )(e- 2aN5 -1), A__ = F(l - We a )(l - e 2aN "), (C.7) 
B++ = e- 2a ( JV5+1 U__ J B— = e 2a ( N * +r >A ++ , (C.8) 

A_+ = ^+_ = 5_+ = £+_ = FW(e a - e~ a ), (C.9) 

C 

~ e aA r s(i _ VFe Q ) - e" aAr s(i _ ITe~ a )' ^ ' ^ 

This notation is valid for positive W and for negative case we define 



e ±a = cosh a ± \/cosh 2 a - 1 (C.ll) 

and flip their sign e — > — e according to sgn(JT). 

The physical quark propagator in momentum space is defined by picking up the bound- 
ary components 

(q(p)q(-p)) = {P L 8 Stl + P R S S , N5 ) ( l ) (8 tyN5 P L + 5 tA P R ) 

V^dwf(P) / 8>t 

= -i ltx sm Pfl G R (N 5 , N 5 ) + WG R (1, N 5 ). (C.12) 

Ignoring the next to leading term in N$ the quark propagator has a simple form 

<9(P)5H0> = ^^- (C13) 
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